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Abstract: In this paper we initiate the study of holographic quantum liquids in 1+1 dimen- 
sions. Since the Landau Fermi liquid theory breaks down in 1+1 dimensions, it is of interest 
to see what holographic methods have to say about similar models. For theories with a gapless 
branch, the Luttinger conjecture states that there is an effective description of the physics in 
terms of a Luttinger liquid which is specified by two parameters. The theory we consider is 
the defect CFT arising due to a probe D3 brane in the AdS Schwarzschild planar black hole 
background. We turn on a fundamental string density on the worldvolume. Unlike higher 
dimensional defects, a persistent dissipationless zero sound mode is found. The thermody- 
namic aspects of these models are considered carefully and certain subtleties with boundary 
terms are explained which are unique to 1+1 dimensions. Spectral functions of bosonic and 
fermionic fluctuations are also considered and quasinormal modes are analysed. A prescrip- 
tion is given to compute spectral functions when there is mixing due to the worldvolume 
gauge field. We comment on the Luttinger conjecture in the light of our findings. 
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1. Introduction 

AdS/CFT [1] is a powerful tool in extracting information about the strongly coupled limit of a 
conformal field theory. The correspondence has been extended to the finite temperature limit 
[2]. In particular, there is recent interest in understanding the thermodynamics, transport 
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Figure 1: In higher dimensions, Landau Fermi hquid theory advocates description using independent 
quasiparticles. In 1+1 dimension, collective excitations are crucial and the quasiparticle description 
breaks down. 

and spectral properties of strongly coupled low dimensional systems, which are of interest in 
condensed matter physics [3-11] (and [12-14] for recent reviews). Quantum liquids in one 
dimension have non-Fermi liquid properties and are thought to be of some relevance to the 
non- Fermi liquid properties of quasi-lD metals [15] and high-Tc superconductors [16]. 

Apart from its ubiquitous appearance and wide applications, one-dimensional Fermi liq- 
uids are of special theoretical interests on their own right. In one dimension, Fermi liquids 
behave very differently from their higher dimensional counterparts, which, in the weak cou- 
pling limit, can be very generally described by the Landau Fermi liquid theory (see for example 
a textbook introduction in [17]). The Landau Fermi liquid theory asserts that the free fermion 
description of the system is only mildly altered under the introduction of interaction, particu- 
larly if we are sufficiently close to the Fermi-surface. Interaction between the fermions can be 
accounted for by introducing a renormalised mass for the fermionic particles. These fermionic 
excitations decay in time, but sufficiently slowly as we get closer to the Fermi-surface, so that 
they are reasonably well defined quasi-particles within their life-time. This picture, however, 
completely breaks down in one dimension. The fundamental difference stems from the differ- 
ent topologies of the Fermi-surfaces, which, in one-dimension, comprises of only two distinct 
points in momentum space, instead of forming a continuous surface. This means that for 
an excitation from the ground state with a sufficiently small momentum, the energy of the 
excitation is completely determined. In higher dimensions a particle hole pair can acquire a 
continuous range of energies even for arbitrarily small total momentum. In ID, this means 
that for small momentum change about the Fermi-points the density fluctuations, which are 
collective fluctuations, are eigenstates of the Hamiltonian. In fact, one can assume a quasi- 
particle description in one dimension and calculate the particle decay rate, which diverges for 
arbitrarily small but finite interactions, indicating a break-down of the quasi-particle picture. 

The Luttinger liquid is introduced to model the effects of interaction in one dimensional 
fermion system, in an analytically controlled context [18]. For a spinless fermion, and ignoring 
back-scattering, which can be shown to be irrelevant in the low energy limit, the model is 
reduced via bosonisation, to a free theory that can be exactly solved, which gives us physical 
insight into the low energy spectrum [19-21] (and for a review see [15,22]). The energy 
eigenstates, as anticipated above, are related to density fluctuations of the fermion liquid by 
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a Bogoliubov transformation. This is a massless excitation with hnear dispersion, which is 
one of the prominent features of the Luttinger hquid. In the more sophisticated model of 
spin-half Luttinger liquid, the diagonalisation of the Hamiltonian leads to two independent 
eigen-bosonic fluctuations, corresponding to the spin and charge densities fluctuations. They 
are again massless excitations which often move with different speeds, depending on the 
coupling. This is the famous spin-charge separation effect [15,21,22], which has recently been 
observed in experiments [23,24]. 

The Luttinger conjecture states [15,21] that any ID model of correlated quantum particles 
(bosons or fermions) having a branch of gapless excitations will have as its stable low-energy 
fixed point the Luttinger model. The asymptotic low energy properties of the degree of 
freedom of this branch will be described by an effective renormalized Luttinger model char- 
acterized by only 2 parameters: a renormalized Fermi velocity v and a renormalized stiffness 
constant K. While successful in providing lots of physical intuition, the Luttinger liquid 
model involves lots of simplification, without which analytic control would have been im- 
possible. In fact lattice methods are used to extract effective parameters of the continuum 
model. The holographic model is thus another analytic handle towards understanding low- 
dimensional systems. The usefulness of these models is that the dynamic properties can be 
extracted more easily than lattice methods. For example, in higher dimensions the fermion 
sign problem makes it next to impossible to extract useful information about fermion corre- 
lation functions using the lattice method. 

The model we will study in this paper is a 1 -|- 1 dimensional defect CFT (dCFT). This 
is made from intersecting D3 branes and is the finite temperature generalization of [25]. The 
counterterms needed to compute correlation functions at zero temperature in this model were 
discussed in [26]. In order to be in a gapless phase, we will consider the phase described by the 
so-called black hole embeddings [27]. Moreover, we will turn on a finite chemical potential so 
that the black hole embeddings are the only physical ones [28] . Considering a defect CFT of 
this kind is also useful since it facilitates comparison with such theories in other dimensions. 
The crucial features of the D3-D7 model with a finite chemical potential which makes a 3-1-1 
dimensional theory with fundamental matter are: 

• A phase transition for low baryon densities and low temperatures [28]. 

• In a certain momenta regime, quasiparticles become a good description for the excita- 
tions in the spectral functions [29-31]. Spectral functions typically contain several Breit- 
Wigner resonances. A collective mode resembling the zero-sound was found [32,33]. All 
modes have a dissipative part. 

• Dispersion relations for quasiparticles look like a; ~ a + bk^ for low k and uj ~ Vck for 
high k. At higher k, the peaks dissolve and the notion of quasiparticles is lost. Vc is less 
than unity and a notion of a limiting velocity emerges [30,34]. 

These features are expected to hold even in the 2+1 dimensional dCFT [35]. 
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The 1+1 dCFT at zero temperature and chemical potential is obtained by placing Nf 
D3 branes in the background geometry created by Nf. D3 branes and considering the probe 
approximation Nf <^ Nc- The probe sees a AdSs x geometry. This background is known to 
preserve half the original = 4, d = 4 supersymmetry and realizes a two dimensional (4, 4) 
supersymmetry algebra. The conformal symmetry group is SL{2, R) x SL{2, R) inherited from 
the AdS^. The massless open string degrees of freedom correspond to a pair of = 4 Yang- 
Mills multiplet coupled to a bifundamental (4, 4) hypermultiplet living at the intersection. 
We send Nc to infinity keeping Nf/Nc ^ 1. This leaves a single M = 4, Yang- Mills multiplet 
coupled to the bifundamental at a 1-1-1 dimensional defect. We will be concerned with the case 
where the probe sits at the origin of the part of the AdS space transverse to its worldvolume. 
The modes corresponding to contracting the inside the saturates the Brietenlohner- 
Preedman bound for scalars in AdS^ and hence are stable. The generating function for the 
field theory is given by the classical action of superstring theory on AdSr, x coupled to a 
Dirac-Born-Infeld theory on AdSz x S^. 

We will be considering this setup at finite temperature. Finite temperature is introduced 
by replacing the AdS^, by the 5-dimensional AdS-Schwarzschild black hole. This trick has 
been used to analyse the behaviour of higher dimensional defect CFTs at finite temperature 
[27-31,35-40]. Part of this exercise is to see if we can probe the validity of the Luttinger 
conjecture using AdS/CFT. The other part is to probe differences between 1-1-1 dimensional 
theories and higher dimensional ones. To introduce a finite chemical potential and quark 
condensate, we look for a probe brane solution in the AdS black hole background with a non- 
trivial brane profile and world- volume electric field. The thermodynamics of the system is 
studied in detail. In the case when the probe wraps the maximal circle (in the notation used in 
this paper, % = 0) we analytically derive the expressions for the thermodynamic quantitites, 
where we find a heat capacity scaling as T in the high temperature limit, as expected in a 1-1-1 
d system, but in the low temperature region where the baryon density becomes important. 
For the non-maximal case (x 7^ 0) we turn to numerics. The specific heat in this case, closely 
resembles that of the X = case. We then compute holographically the Green's functions of 
various scalar, vector and fcrmionic operators by studying the world- volume fluctuations of 
the probe brane about the background solution. 

Our analysis is divided into several steps. To begin with, we study the fluctuations of 
the longitudinal electric field for trivial embedding (where the brane passes straight through 
the horizon), very much in the spirit of [32]. We find that there is a massless (zero-sound) 
mode with dispersion given by the conformal result i.e. uj = ^k/^/q, where q is the defect 
dimension, equals to one in this case. We obtain this result analytically in the hydrodynamic 
limit, and find no dissipation. We extend the study to larger frequencies oj and momenta k 
numerically, and find that the dispersion relation is unmodified, and remain dissipationless. 
More surprisingly, while it is by now known that these massless (sound) modes in holographic 
defect models disappear at finite temperatures for higher dimensional defects [41], the massless 
mode in our ID defect survives at finite temperatures, with identical dispersion as in the case 
of zero temperature. 
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We then generalise the investigation to the case of mixed excitations, where the fluctu- 
ations of the electric field longitudinal to the defect are mixed with those of the embedding 
profile at finite momentum exchange, and in the presence of background world- volume electric 
field and a non-trivial probe embedding. To deal with the mixing, we explain in detail how 
the quasi-normal modes and spectral functions are obtained in principle [42], and implement 
the procedure numerically. It is interesting to see that the massless mode mentioned above 
survives even in this limit without dissipation. There is also a mode corresponding to pure 
dissipation, which originates from the profile fluctuation even before mixing is introduced. 
The dispersion of this mode is given schematically by a; ~ —i{a + bk'^), for some constant 
a and b dependent on the background electric field (controlling the chemical potential) and 
the probe embedding (controlling the quark mass and condensate). The higher quasi-normal 
modes for bosonic excitations, however, behave roughly as the higher dimensional ones, with 
a dispersion relation quadratic at small lo and k, and approaches that of the speed of light 
asymptotically for large uj and k. 

It is interesting to note that while our analytic expressions for the dispersions for the 
massless mode appear to be independent of the charge density d, the limit d — > at zero 
temperature is not a smooth one. In fact when d = the fluctuations can be solved exactly 
analytically at zero temperature and we flnd that the massless mode disappears. In fact the 
spectral function becomes a constant, for both the vector and scalar modes. At zero d and 
finite temperatures, the scalar fluctuations can again be solved exactly and the quasi-normal 
modes form a discrete infinite tower of modes. The massless dissipationless mode that appears 
in the longitudinal electric field fluctuations however, remains intact here. 

This paper is organized as follows. In section 2, we discuss the setup and specify our 
conventions. In section 3, we consider the thermodynamics of our embeddings and carefully 
analyse the counterterms needed for the calculations. In section 4, we turn to the analysis of 
zero sound. In sections 5 and 6 we consider the spectral functions for bosonic and fermionic 
fluctuations. We conclude with a discussion in section 6. Appendix A gives some analytic 
results for the thermodynamics. Appendix B has a proof that when the real part of the 
quasinormal mode is non-zero, then the imaginary part of the mode has to be necessarily 
positive indicating no instabilities. Appendix C has some explicit calculations for the Green 
functions for the mixed modes considered in section 5. 

2. D-brane configurations 

We would like to begin with a review of the brane configuration, namely, the intersecting 
D3 systems, considered in this note. Its zero temperature limit has been studied in depth 
in [25,26]. We have Nc D3 branes and Nf D3 branes intersecting over a 1-1-1 dimensional 
domain with four relatively transverse dimensions. The low energy effective field theory is 

given by a supersymmetric U{Nc) x U{Nf) Yang-Mills theory on a 1-|-1 dimensional defect. 
In the large Nc and large t'Hooft coupling A = At^QsNc limit, such that Nf/N^ <C 1, one 
could replace the Nc branes by the AdS geometry and treat the Nf D3 branes as probes 
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in the curved background. The excitations on the probe are then dual to the low energy 
excitations on the 1+1 dimensional defect. Although the 1+1 dimensional defect will be the 
primary focus of this paper, sometimes we will draw parallels to the well studied D3-D7, 3+1 
dimensional setup as well. The brane-scans for the two setups are 
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NcBS: xxxx 

1. Nf D7: xxxxxxxx 

2. Nf B3: X X x x 

In the zero temperature limit both these theories are supersymmctric. We will consider 
the abelian case only and leave the analysis for non-abelian effects for future work. At finite 
temperatures, the AdS geometry is deformed to that of an asymptotically AdS black hole [43] . 
The metric is given by 

3 
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where the metric for the 5-sphcrc can be chosen to be 

dnl = d9^ + cos^ edCf + sin^ e{dCi + sin^ (2dCi + cos^ (2(^(4), 



(2.3) 



The temperature T of the system is related to the radius of the black hole horizon rn by 



T = 



7rL2 



(2.4) 



Embedding 

In the D3-D7 setup, the probe wraps {xq, xi, X2, X3, u, ^2, Csj C4} while in the D3-D3 setup 
the probe wraps around {xo,Xi,u,(i}. We will turn on a world- volume gauge field Ao{u) 
which corresponds on the gauge theory side to having a chemical potential or a finite baryon 
density. In the presence of a chemical potential, it can be shown that the so-called black 
hole embedding is a physical embedding, which we will focus on. It is common to make a 
change of coordinates cos 6 = x(«) in the case of D3-D7, and for D3-D3 it is convenient to 
pick instead sinO = x(^)- For black-hole embeddings, we require regularity at the horizon. 
In the literature [30], this is implemented by switching to a coordinate p = x/u^ + V'w^ — 1 
and requiring that dpx{p = 1) = O"*^. We also set x('u = 1) = xo- 

The resultant probe brane action for the D3-Dp setup is given by 

^ = 4uik = J df+VVdet(G + 27r«'F) 
= J d'^aduL^ovik 



dVd?x(l - X 



p-i 
4 
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(2.5) 



^In the near horizon hmit p ~ 1 + \/u ~ 1. Therefore requiring dpx{p = 1) = is equivalent to requiring 
that the expansion of x in m to go like xo + Xi(w — 1) + 
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where g is the induced world-volume metric and explicitly for the D3-D3 case, and allowing 
for a non-trivial profile in xi^), we have 



L2 



-{I - \)dt^ + dx'^ 



2 



{l-^)-\l + u'iP^)du' + u\l-x')dQl . 

«4 (1 - X^) 

(2.6) 

Here, F = dA are world-volume gauge field strengths. In the Euclidean signature, the integral 
over the time direction is between zero and 1/(^^7") = 1/T in our units where ks, the 
Boltzmann constant, is set to one. This factor is made explicit after the second equality in 
equation ( |2.5| ). The normalization constant N is given by 

where q is the spatial dimension of the intersection domain, and Vp-q-i is the volume of the 
sphere wrapped by the probe Dp brane. We choose the gauge Au = 0. The field strengths 
can be readily solved in terms of x- 



guu\m\{(i^ - l^ool^ii) 



9ii{b'^9ii - (P\9m\ - (511(1 - x^)) Iffool) 



(2. 



where 6, d and e are dimensionful integration constants from solving the equations of motion 
of 

Fux^Fqu and Fqx respectively. They in fact parametrise the conserved magnetic and electric 
charges of the solution. Note that we have absorbed a l/27ra' in the definition of the constants 
6, d and e. For notational simplicity we will take L = 1, unless stated otherwise. 

We will concentrate on the simple case where there is a finite baryon density, correspond- 
ing to a non-trivial d, while leaving b = e = 0. 

f id/r%) _ 

Ao{u) ~ <^ ^+ 9- (2.9) 



-|- d log u q = 1. 

near the boundary. The thermodynamics interpretation of various quantities will be discussed 
in the next section. So for the numerics we have two tunable parameters xo a-iid d. The near 
horizon expansion for x is given by 

X = X0 + Xiin-l) + O{{u-lf), , (2.10) 

with XI = ^5°^^ ~ for D3-D7 and xi = ^^"(^ ~ ^0) for the D3-D3 setup. 

4(d2 + (l-xg)3 4(d2 + (l-xg)) 

Here we have introduced dimensionless d defined as 

d=-ir- q = 3, d=— q = l. (2.11) 
Typical profiles for D3-D3 at various values of d are plotted in figure (§). The asymptotic 
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expansion for x reads 

m c 

^logn (2-12) 

\ 9=1- 

u u 

In the supersymmetric limit m/u is an exact solution, both for q = 1 and q = 3. Hence, we 
will continue to use the same notation m and c for the D3 probe, as in the higher dimensional 
case. At q = 1, X saturates the Brietenlohner-Freedman bound. The asymptotic behaviour 
presented in ( p. 121 ) is as expected since the conformal dimension of the dual operator A± = 
(q + l)/2 ± y^(g + 1)2/4 + m2, where 771 is the mass of the scalar field {fh? = — 1 for x): 
evaluates to A_|_ = A_ = 1 for x at g = 1. This already suggests that the two possible 
large r behaviour are log(r/ro)/r and 1/r, as we have confirmed. This however, also implies 
that only one conformally-invariant boundary condition is possible [44]. One could require 
physical fluctuations to behave as 1/r and setting the log(r/ro)/r terms to zero. These, as 
mentioned above, correspond to the supersymmetric solutions. On the contrary, requiring 
the vanishing of the 1/r term at the boundary and keeping the log(r/ro)/r term violates 
conformal invariance, since that is a scale dependent statement i.e. requires the choice of an 
arbitrary scale rg. There is therefore likewise only one natural way to incorporate an external 
source cq, and that is to require that x behaves like cq log(r/ro)/r in the AdS boundary 
r oo. From the point of view of the dual field theory, this means that the expectation 
value of the dual operator has a logarithmic violation of conformal invariance in the presence 
of a source, since the two point function contains a logarithmic divergence. Considerations 
here apply similarly to the leading log term of the gauge field in (^]^). We will encounter this 
issue again in the discussion of the thermodynamics of the system. 

Therefore, the roles played by m and c are interchanged, namely that m corresponds to 
the quark condensate and c to a source [26]. This is in contrast with the q > 1 scenarios [27]. 
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We will make further comments about the appearance of the log as we proceed. 



3. Thermodynamics 

In this section we consider the thermodynamics of the 1+1 defect CFT. There are some 
interesting differences as compared to higher dimensional defects as will become clear. These 
differences arise due to the logtt terms in the asymptotics of x ^'^'^ ^t- We introduce the 
dimensionful source C = ckT, which is the ID analog of the quark mass in higher dimensional 
probes [27], and 

, - j dCi^^ = 27r{27Ta' NfTD3)d. (3.1) 



n, 



Since C and Ug are parameters of the dual defect CFT at q = 1, the phase diagram of the 
system is obtained by keeping them fixed as the tcmpcratmc T is varied. The ratio Ug/C is 
therefore a handy parameter to label the dual CFT. We will find that below a certain value 
for Uq/C the system exhibits an instability, similar to the D3-D7 case. We will be interested 
in values of Uq/C where the system is stable. It should also be noted that while the Mermin- 
Wagner theorem states that at finite temperatures continuous symmetry cannot be generally 
broken in dimensions d <2, the restriction is evaded in the large N limit [14]. 

To analyse the thermodynamics of this system we begin by computing the Euclidean 
action on-slicll. From the asymptotic behaviour of the fields, it is clear that the action is 
divergent. The counterterms needed to cancel the divergence arising from x have been given 
in [26]. 

Li = -y Vt, (3.2) 



4 



L2 = -log(A/ro)V7i?^, (3.3) 



L4=yV7X^(x,A)(l-j^^^), (3.4) 

where 7 is the induced metric on the boundary, is the Ricci scalar evaluated on 7, and 
we have defined A = A/rn- Since we are considering a flat boundary theory, L2 does not 
contribute. It is important to note here that, since we have switched coordinates to dimen- 
sionless u, all functions of u in the boundary limit is evaluated at A. On the other hand, 
when the cut-off appears explicitly as a coefficient in the counter terms, it appears simply as 
A, and in the case of explicit log terms, log(A/ro), for some arbitrary scale rg. To make sense 
of the dependence on ro, one should recall that there is a logarithmic violation of conformal 
invariance at non- vanishing d and c, as discussed in the previous section. Therefore there 
is dependence on an arbitrary scale ro, which appears inside these log's. The dependence on 
logro reminds us that whenever we perform a rescaling in r i.e. r Ir, corresponding to a 
rescaling to a different energy scale in the dual theory, a contact term (which contributes only 
to the finite part of the Green's function) proportional to log I would appear in the action, as 
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discussed in [45]. This will not affect the physics (e.g. spectral functions, quasi-frequencies 
etc) we are interested in and we will set rg to unity in the rest of our discussions. 

For non-vanishing d, there is also a logarithmic divergence in the action of the form 
-d721ogA. This is a new divergence only appearing in 1+1 dimensions. We will add 

^F = ^^.A.7^'^^, (3.5) 

as a counterterm to remove the logarithmic divergence due to the gauge field. Note that the 
term preserves gauge invariance despite looking otherwise. The reason is that AdS/CFT puts 
restriction on the allowed gauge transformation 6A = d$ such that the asymptotics of 
are not altered. Since A^|„^oo ~ + ^ogu, the leading term at the boundary of could 



at best be n . The gauge transformation of (3.5) again goes like u and thus vanishes at 



the AdS boundary. It is shown in the appendix that in the gauge Au = 0, the boundary term 
can again be written in terms of gauge invariant variables. 

What is the thermodynamics interpretation of n, d, m and c? Firstly we observe that with 
a radial electric field, there is effectively a number Uq of fundamental strings stretching along 
the worldvolume of the probe D3 brane. This density of strings is given hy d = —5Iy,uik/SFtu- 
As a result it is natural to interpret d as being proportional to a number density. This 
leads to identifying fi as the chemical potential^ . Similarly in cases of the higher dimensional 
defects, the brane separation which is controlled by m is identified as being proportional to 
the inverse temperature. However, unlike the higher dimensional defects, for the probe D3, 
m could be interpreted as the brane separation only in the supersymmetric theory. Moreover 
from the point of view of AdS/CFT, the source of the dual operator should be identified as 
the coefficient of the logu/u term i.e. c. As a result we will expect 

This is consistent with the interpretation in [26] where m is identified as a vev in the dual field 
theory. Thermodynamically, we then have that the Euclidean action is a function of Uq and 
the source c rather than [i and m as in the higher dimensional defect case. Since C = c{tiT) 
is kept fixed and corresponds to a source in the dual defect CFT, we can identify 1/c as the 
temperature. The free energy T can thus be interpreted as the Helmholtz free energy. 
The counterterms arising from Li + L2 + + L4 computes to 



Let 



d'^a[ - — H ^ + -(m + clogA)"' - mc 

^ ^ ~ {d^ - r'jjCpyLogru 

and we define T jk^T = I = /buik + ^ct- Note that the last term above arises precisely because 
of the explicit appearance of log A/ro in the definition of the counter terms. 

^As we will see below the physical chemical potential Ji differ from /i by a constant in the case of a probe 
D3 brane. This applies also to the discussion of m, which will be discussed below. 
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After considering the variation of the action, this leads to 

/N = —{m — c log rnJSc + {/J. — d log rn) / rnSd. (3.8) 

This confirms our expectation that the action is a function of c, d. It is not surprising that 
m — clogrn and /x — dlogrj^ appears in this particular combination. This is because the 
dimensionful vev's are defined as 

X(r) = -CM +Clogr), Ar(r) =JI + dlogr, (3.9) 
r 

implying that M + Clogrn = rnm and similarly /I + dlogr h = ^. It is thus consistent to 
have these physical quantities M and /I appearing in the variation of the free energy. 

3.1 Specific heat 

Having discussed the interpretation of the Euclidean action in the previous section, we are 
ready to compute the specific heat capacity. 

3.1.1 Case I: x = 

For simplicity, let us begin with the massless case where x = 0- The entropy is obtained by 
diflFerentiating the action at constant chemical potential. The chemical potential is given by 
the constant in the asymptotic expansion of the gauge potential ^o- 

Ji = II — dlogrn = lim Aq(u) — AQ(u)(ulogu) — dlogr h, ^oiu) = / duF^Q. (3.10) 

At X = 0, they are explicitly, 

F„o = A'oiu) = ''"'^ , JI = d(log 2 - logirn + \J(P + rl)). (3.11) 
J rjjU^ + 

Combining with the counter terms, we have 

= ^{-2rH^d^ + rjj + d\l + log 4) - 2d^ log{rH + sjd? + r]j)). (3.12) 

As discussed in the previous section, at zero Xi ^ = ksTI is simply a function of rn and the 
baryon density d. Therefore the entropy S is given by 

S_ l_dT 

iVo ~ ~l^dT 

and the heat capacity Cy is 



TT^/d^ + (7rr)2, (3.13) 



c^^^d{S/No) 



No dT 



<'^'^ (3.14) 



J V(tT)2 + <P 



The heat capacity is linear in the temperature at suflficiently high temperatures as expected 
of a 1 + 1 dimensional quantum system. However, the temperature dependence becomes 
quadratic at low temperatures in the presence of a non- vanishing d. Note also that the result 
is unchanged if we subtract the zero temperature contribution to the free energy and the 
chemical potential. Note that in this limit i.e. x = 0) the entropy and subsequently the heat 
capacity are independent of the arbitrary scale ro, if we are keeping the baryon density fixed. 
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3.1.2 Case II: x / 

When X is non-zero, the above calculations can be done numerically. However, we would then 
have to take into account the contribution of x s-^id it's corresponding counter terms. We 
follow the procedure in [27] to obtain an explicit expression for the entropy S = —ndr^jJ-H, 
keeping the baryon density d and C = rnc fixed. i.e. 

c - d 
dr„c = , drijd = . (3.15) 

The entropy thus has five contributions. 

S = -'K{Si + Sii + Siii + Siv + Sv), (3.16) 

where 



F -A 
Si = (a,.^AA) — , Sii = ( )Lbuik{u) 



A' 



S, = -M ^ , ' dr„ logTH = -M '\^ ' ■ (3.17) 



Siii = ( — —)dxlct, Sijj = (— )(c9c + dd/)J^ = —^{dn/rH - cm + (cP - c^) log rn) 
rn rn rn 

We have made used of the fact that jN = — (m — clogri^)(5c + (/u/r// — d\ogrH)&d in 
evaluating Si^, and S^ is the explicit derivative with respect to ru in the last {d^ — (?)/2 log ru 
term in the counter terms in equation ( |3.7| ). 

After a tedious but straight forward computation, we finally have 



'5 = -2^(1 



[fJ'/^H — dlog{'irT))d — c{m — clog(7rT)) 



2T ' 2T 

The corresponding heat capacity is given by 



^J: ' . (3.18) 



TdrS 



= S + ^dTT\{fi/rH - d\og^T)d- {m - clogTTT)c) - M^-^--^. (3.19) 



Derivative of c and /i with respect to T has to be implemented numerically. Typical plots 
of the free energy against temperature at fixed d/C in the stable and unstable regimes are 
shown in figure (|3|)^. At high temperatures, both the entropy and the specific heat scales as 
T, approaching the conformal result. Below d/C <1, the heat capacity could turn negative. 
This allows us to determine the phase diagram of the system. This agrees with a quasi-normal 
mode analysis, where an unstable mode that grows in time appears where the heat capacity 
is negative. The phase diagram is shown in figure (Q)^. 

^The dotted line is added in by hand. Our numerical results end just before the dotted line begins and the 
region corresponds to very small d ~ 0.01 and an initial condition for x at the horizon extremely close to one 
Xo ~ 1 — 10~^. It would require much more numerical accuracy to see the dotted line. 

''At T = 0, d = 0, the supersymmetric solution x = n'*/^ is stable. Note that one could not smoothly go 
from a general black hole embedding to a supersymmetric solution by taking the c — > limit. For more general 
X,the Minkowski embedding is the only stable embedding [27]. 
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Figure 3: (a) and (c) are plots of free energy and the corresponding entropy against temperature in 
the stable regime, where d/C > 1. (b) and (d) are typical plots of T and S in the unstable regime, 
where d/C < 1. The dotted line is added in by hand. 



4. Sound modes 



For simplicity, we would like to begin by considering the simpler embedding where x = 

all the way. The terms in the action mixing the gauge potential and the brane profile 
perturbations are zero in this case. Therefore we can consider the equations of motion of 
the longitudinal gauge perturbations independently of the other scalar fluctuations. The 
calculation in the zero temperature limit mimics that in [32]. The solution there however 
cannot be applied directly for probe D3's by simply substituting the defect dimension q = 

1 in our case. The probe D3 case however is straightforward to solve and one obtains a 
dissipationless mode with a dispersion given by the conformal value to = zizk. In the finite 
temperature, the near horizon solution has a singularity and requires slightly more work. It 
is convenient to work with the coordinate z = 1/u, which has the correct normalisation at 
the boundary and would avoid mistake when we begin to expand in power series of to and q. 
The equation of motion at zero x is given by 
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Figure 4: The phase diagram of the system. The shaded region (excluding the origin) is unstable. 



d,iF[z]E'^[z]) + G[z]E,[z]=0, 

^ z(l + d^z^)3/2(-l + z^)u; 

2g2(-l + z4)+ 2(1+^2^2)^2^ 



z' 



G[z] = ' — (4.1) 

^ ' 2(-l + z4) ^ ^ 

We would like to solve these equations in the limit, where uj and momentum k are much 
less than any other scales, in which case we can obtain a solution perturbatively in to and k. 
We will follow the approach in [46]. To begin with, we isolate the singularity appearing in 
the near horizon limit, and write 

E^[z] = {1- z)^{eo[z]+ujei[z] + ...) (4.2) 



After making the substitution in equation ( |4.1| ) and removing the overall factor of (1 — 
z)^^ , we expand the resulting equations of motion in powers of to and k. In the expansion 
we considered the case u; ~ fc. For each solution of eo,ei... we have to set two boundary 
conditions. Since the singular part of the solution has been isolated, we require that eo[z] 
and ei[z] and so on, are regular at the horizon. Also, by convention we choose to absorb any 
constants ei>o[-z 1] by eo, and therefore require that ej>o[^;] vanish at the horizon. Given 
these conditions, we find that 

eo[z] = Co 

ei\z\ = co(ci 5 h other terms) (4-3) 
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The solution ei{z\ involves other complicated terms in log[l + (Pz^] and log[l + d?z'^ + 



(1 + d?){l + (Pz'^y\, which are regular at the boundary and do not concern us here. The 
constant c\ can be expressed in terms of d, a; and fe, by requiring that ei vanishes at the 
horizon. What is important however, is that the coefficient of the Iog[2;] term is given by 
(fc^ — o;^), which gives a quasi-normal mode at 

uj = ±k, (4.4) 

in agreement with the conformal result of a sound mode in one-dimension. We can in fact 
evaluate the spectral function explicitly. To lowest order in lo and k, it is given by 

lm{Go;x) Co Lv^ 

rM— -TTCTV) = T-9. (4-5) 



(2Ar)/(7r2r2) a;2-fc2 ^ icojk^-w^) ^ (a;2-fc2)2' 

and similarly for other components of the Green's functions of the currents, which are obtained 
by putting into the above expression appropriate pre-factors of /c's and cj's. 

The real part of Gxx approaches a finite value in the hydrodynamic limit at = 0. This 
is because as we consider higher order terms we find 

hm Re(G..) ~ Re(-^3^^^^^^-^^^^^) 

= Hi-Ki (4.6) 

where Hr i and Knj are overall zeroth order constants depending on w, k. We will have to 
solve for still higher order corrections in order to determine Ki. 

We followed this sound mode using the complete equation without taking the hydrody- 
namic limit and found that this dispersion relation remains exact beyond the hydrodynamic 
limit. In fact, even when x is non-zero and mixing becomes important, as we will investigate 
in the following sections, this mode remains intact and dissipationlcss. 

This should be put in sharp contrast with higher dimensional defect theories obtained 
from higher dimensional probes (i.e. D5 or D7), where the leading term in the boundary 
expansion of the gauge field is given by the constant term. However, for regularity at the 
horizon, eo continues to be a constant in these higher dimensional theories. The constant 
term at the boundary is then given by 

Ex[z]=co{l+oo{ci) + ...) (4.7) 

where cq is an arbitrary constant, and ci depends on k and a; although it is overall zeroth 
order in uj and k. The quasi-normal modes are defined at the zeroes of the constant term. It 
is however clear that uj cannot begin at linear order in k, due to the presence of the constant 
term cq. This absence of a sound mode has been observed in [41], and is interpreted as the 
destruction of the Fermi-surface at finite temperature in the strongly coupled regime. 

We have also looked for quasi-frequencies of the Q fluctuations in the hydrodynamic limit, 
and for simplicity, at zero x ^-^d zero chemical potential. In that case the regular solution 
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zeroth order in to is non-trivial and contains a log term at the boundary, whose coefficient 
is independent of u. The quasi-normal mode occurs at the zeroes of the coefficient of the 
log term, and therefore no higher order corrections in uj and k could cancel the zeroth order 
contribution, if we assume that to ^ k. We conclude that it does not contain a sound mode 
with linear dispersion. 

While the above expression for the quasinormal frequency appears to be independent of d 
that controls the chemical potential, the sound mode actually disappears at zero temperature 
if we then also put d = 0. The limit to d = is thus not a smooth one. The difference lies 
in the modified near horizon (i.e. u ^ or z — > oo) expansion of the solution. Instead of 
exp{±iijz)/z we have ex.p(±iz\/uj'^ — k'^^j^/z. We pick the plus sign for regularity as in [32]. 
The equations of motion at d = and T = can in fact be solved exactly. The solution is 
given by 

^^(a;, z) = ciIo{z^/k^-uj'^) + ^K^{z^/W^), (4.8) 

where / and K are the modified Bessel functions. The constants ci and C2 are fixed by the 
boundary condition, and they are related by 




(4.9) 



Then we can expand the solution in the boundary limit, and extract the coefficient of the log 
term and the constant. They are given by 



£;^(z^0) = i7 + 7r + ilog(^-^ ) +ilog(z) + ... (4.10) 

where 7 is the Euler's constant. The log term becomes independent of lo and k and the 
massless sound mode disappears. 

We can also obtain the Green's function from the solution^. At zero k this simply gives 

G.. = ^(-7 + ^vr-log(|^)). (4.11) 

A plot of the real and imaginary parts of the Green's function are shown in figure (^. The 
real part depends logorithmically on w, where as the imaginary part is simply a constant 7r/2. 

The result also agrees precisely with the high lo limit of the numerical solution at finite 
baryon density and temperature, which we will return to later along with the numerical 
solutions. 

^We obtained the results at th ~ 1- As we will show more explicitly in the appendix, the explicit log rj/ 
contribution in the counter term contributes only to the real part of the Green's function and does not affect 
the imaginary part, which is of more interest for our purpose. 
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Figure 5: The red and blue curves plot respectively the real and imaginary parts of G^x at zero d, T 
and k. The imaginary part is a constant given by 7r/2. 

5. Bosonic fluctuations and spectral functions 

5.1 The diagonal modes 

Consider the fluctuations^ of the angular position C,2^ , as defined in section (^) of the probe 
D3 brane. The diagonal modes exhibit the usual quasi-particle peaks and indicate no sign of 
the underlying instability in the system. Here we plot some representative spectral functions 
for the D3-D3 setup in figure (P), since it has never been studied before. The D3-D7 case was 
studied for example in [30]. The first peak is located precisely aX uj = k. The width of the 
peak increases with A;, and the Green's function has a singularity at precisely w = /c = 0. 

5.2 The mixed modes 

In the presence of a chemical potential and at a finite momentum, there is a mixing between 
the fluctuations of the longitudinal electric field Ex and the profile 56. This was first pointed 
out in [30]. These mixed modes were subsequently studied in detail in [31] for the case of 
probe D7's. This Lagrangian at quadratic order can be written as 



+ Di6eA + D25euFtu + D^SOFtu 
+ S^Fl + S 2Fl + S^F^xu- (5.1) 

"Note that we have not considered fluctuations of the probe along x^,x^ inside the AdS^. The action of 
these modes include contributions from the WZ terms [25,47,48], and the boundary action evaluates to zero. 
They are related to surface operators in the dual theory [47,48]. The usual mesonic excitations in the dCFT 
are described by fluctuations of the gauge fleld and of the transverse 3-sphere [27]. We will not discuss these 
surface operators any further in this work. 

'^We have not considered fluctuations of other angular directions on the transverse 3-sphere, since they are 
related to the ("2 mode by 50(4) rotations. 
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(a) (b) 




Figure 6: The real (red curve) and imaginary (blue curve)parts of the Green's functions corresponding 
to C,2 fluctuations, at rn — 1, d/C — 1.074 and (a)fc = 0, and (b) at fc = 2. 



To avoid cumbersome notations, we have implicitly made a change of coordinates i.e. (t, x) — > 
■^{t,x) . As a result our 2-momenta are dimensionless quantitites normalized by the tem- 
perature. Here Si,Ti,Di are complicated functions of u whose complete expressions we will 
refrain from showing here. Their near horizon and boundary behaviour will be given in the 
appendix. The second line denotes mixing terms and all the mixing terms are proportional 
to d. We will work in the gauge = 0. This leads to the constraint equation 

dtD26eu + dtDsSe = 2dt{S2duAt) + 2d,{S3duA,) . (5.2) 

The equation of motion for A^. and Af are 

dt{Di6e.,) = -2dt{SiF.,t) + 2du{S3duA,) , (5.3) 

-d^{Di6e^) + du{D26eu) + duiD^Se) = 2d,{SiR,t) + 2du{S2duAt) . (5.4) 



Now we write 



A,= /^e-*+^'=^a,(n,K), (5.5) 
Se= I ^e--'*+-'-e(tx,K). (5.6) 



Since Afj_ and 59 are real, we have 

= e(u, -k) = 6(n, k) = e^, a^^^^ = af,{u,-K) = afj,{u,K) = a^^^ , (5.7) 

where ~ denotes complex conjugation. The gauge invariant combination is Ex = kat + ojax- 
Using the constraint equation we find 

_ kl^ + 2iuS2E', 
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or equivalently 
where 

A = -iujD2@' - iuDsQ . (5.10) 



Using this we find 



z)..*e = -25,.£. + 2a.|S3|^^^). (5.11) 



After some algebra, the equation of motion for 6 is found to be 

(5.12) 

Note that our equations of motion exphcitly violates spatial and time reversal. This 
is a result of the presence of a background worldvolume electric field F^t- This will have 
implications in the spectral functions, which we will discuss in the next sub-section. 

5.2.1 The spectral function of the mixed modes 

Recall that for the mixed fluctuations, we have two coupled linear second-order differential 
equations, and as a result, there should be four independent pair of solutions (_E,* (n), 0*(n)), i G 
{1, ...4}. By imposing the infalling boundary conditions at the horizon for both Ex and 0, 
we are imposing two constraints and should be left with two independent solutions. The 
boundary values of and respectively source different operators, and given these two 
independent solutions, it is in principle possible to construct solutions with independent ar- 
bitrary boundary values E^ipo) and G(oo). i.e. in general 



Ex{u) = AEi{u) + BE2{u), Q{u) = AQi{u) + BG2(n), (5.13) 

which allow us to solve for A and B for any given pair of boundary values {^^^(oo), B(oo)}. 
The Green's function is obtained by differentiating the on-shell action evaluated at the bound- 
ary with respect to Exipo) and 6(oo). It is therefore apparent that while explicit mixing 
terms in the on-shell action vanishes at the boundary, mixing terms involving the product of 
Q{qo)Ex{co) could still appear through the squares and products of A and B, which generally 
depend on both of these boundary quantities. It is important to emphasise that the sources 
{Ex{oo),Q{co)} are only independent if we have the freedom to pick A and B accordingly. 
Given a particular boundary condition at the horizon, the boundary values of the pair solu- 
tion {Ex{z)^@{z)) of the equations of motion are not independent and therefore evaluating 
the action at such a particular solution and differentiating with respect to these correlated 
boundary values would not give the correct Green's function, as in [31]. In fact, it is apparent 
from the procedure taken in [31] that the diagonal elements of the spectral function (i.e. imag- 
inary parts of the Green's function) are not proportional to the conserved current discussed 
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in section ( ^.2.2 ). By picking different horizon boundary values these spectral functions could 
turn negative, which already suggests that the procedure is pathological. 

The study of mixed modes is well known in the context of shear modes in R-charged 
backgrounds [49]. The numerical procedure needed to search for quasi-normal modes have 
been discussed in [50]. We will discuss this procedure, and lay out the explicit way of comput- 
ing spectral functions^. To implement the procedure numerically, where we are only capable 
of controlling the horizon values, we construct general solutions of and Q in the following 
manner. First, we construct two independent sets of solutions by picking two sets of horizon 
boundary conditions, for concreteness, say 

Ei{u ^ 1) ~ (1 - u)^{l + ...) Gi(n ^ 1) ~ (1 - u)^{l + ...), 
S2(n^ 1) ~ + G2(u^ 1) ~ (l-u)^(l + ...), (5.14) 

where the ellipses again denote higher order terms in (u — 1), which are solved in terms of 
the lower terms. Then we define two pairs of solutions 

ex = Ei{u,k) - ^^^^°°''^I e2{u,k), 60 = Gi(u,k) - ^4-^^^02(n, k) , (5.15) 

+ TP f \ -^1(^00, Z^) p . ^ , ^ f X Ei{Uoo,k) 

tx = El{U, K) - —- -E2{U,K), te = Ol[U, K) - —- -k)2{U,K), (5.16) 

E2{Uoo,K) E2{U^,K) 

where we denote k = {uj,k). These solutions are constructed such that the leading boundary 
terms in eg and vanish. As a result, a general solution can be written as 

Eo{k) f ' \ +Tq{k)^ f ' > ] , (5.17) 

Q = ne[ Eq[k) + Tq{k)— , (5.18) 

where for D3-D3^ 

ne = ln«oo, ng = , (5.19) 

the normalizations chosen to respect the asymptotic boundary conditions. Eq and Tq de- 
termine the boundary values of E^ and Q respectively, and source the corresponding dual 
operators. The boundary action can be shown to be schematically 

-M I dujdk^^^ExE' + TiQQ' , (5.20) 

J U!'' — k"^ 

with which we can compute a matrix of Green functions: 

/ 2 e'^{l^,u) t'^(^^u) ^ e'g{-K,u) 



*We thank Andrei Starinets and Sean HartnoU for discussion of the issue and pointing us to useful references, 
^for D3-D7 ne = l,ne ^ — . 
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where we have abosrbed the factors of — k'^ inside Cx to avoid clumsy notations. We 
also need to include the contribution due to the counter terms^'^. When the divergences are 
removed the Green's function would be given by 



G = M 



,(1), 



(k) 



-k) 



\ 



if\k) 



W(-k) 



tf\k) 



(5.22) 



where we denote with superscript (0) for the coefficients of the leading boundary term and (1) 
for the coefficients of the first sub-leading boundary term in the boundary expansion of the 
solutions. Clearly, the vanishing of the leading boundary term, corresponding to the defining 
condition of the quasi-normal modes, again coincides with singularities of the Green's function 
matrix. It is important to note that in ( ^.16 ), by construction , the zeroes of e^*^^ and t^^^ 
coincide, and are proportional to 



QNM = ^162 - -^261 



(5.23) 



This is the same definition for quasi-normal modes as in [50]. This also implies that all 
the entries in the Green's function matrix have a singularity at the same time. This is well 
known from the study of the coupled shear modes from the fiuctuations of the background 
metric [49]. We will work out the Green's function matrix explicitly for the case of D3-probes 
in the appendix. 

5.2.2 Symmetries of the action and the spectral function 

In the absence of mixing, we have 



MA 



to be symmetries of the action. This leads to a U{1) x U (1) symmetry for each k in the bulk 
action which leads to two (infinite sets of)conserved currents Je^ Je,k- (This is probably 
an accidental symmetry of the quadratic action, to be broken by higher terms.) 

In the presence of mixing, only (pi = (j)2 leads to a symmetry of the action. This happens 
when one turns on a non-zero d and a non-zero k. In this circumstance, the global symmetry 
breaks to a U{1) with a single conserved current. This works out to be 



J{u) = 2Im 



(5.24) 

Note that terms containing are not present in the conserved current since it is ~ |0k|^ 
which is real and is canceled out. Thus, we have a conservation equation 



duJ{u) = 0, 



(5.25) 



■'Again we are setting th ~ ^■ 
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which leads to a condition for checking our numerics. As noted in [46], the imaginary part 
of the boundary terms of the Lagrangian evaluated on shell is proportional to the conserved 



current ( 5.24 ).i.e. 

^J(u ^ oo) 

and is thus independent of the radial coordinate u. Also we should note that terms containing 
Di's (i.e. all the mixed terms) in fact asymptote to zero at the boundary. 
Explicitly, the current evaluated on-shell at the horizon, is given by 



J(l) = Jl + d^-xliu; + -), (5.26) 
where we have made use of the near horizon solutions of and Ftx that satisfy the infalling 



boundary conditions as given in ( 5.14 ). The result is independent of the choice of signs in 

O)- 

The first term in (|5.26| ) originates from the Q@' term while the second term from clxCl'^ 
and ata[. The cross terms again vanish at the horizon. The sign of the expression is explicitly 
positive definite, and this follows from the signs of T^, 5*2 and 5*3, and also the infalling 
boundary condition. In the absence of coupling, U{1) x U{1) symmetry is restored, such that 
the two terms in J(l) are conserved independently. 

Recall that the spectral function are obtained from the quadratic boundary action in 
eq. ( ^.26| ). In fact, from our discussion in ( |5.22| ), the imaginary part of each of the diagonal 
elements of the of the Green's function matrix, is simply given by 

iMGee) = Im(G,,) = (5.27) 

l^ol \To\ 

where J{ex, ej) denotes that the current is evaluated on the first pair of solution {ex, et), and 
similarly for J{tx,tt). Note that the counter terms, which are explicitly real, do not alter the 
value of the imaginary parts of the Green's function. The imaginary parts of the diagonal 
elements are thus proportional to the conserved current, which from ( |5.26| ), is shown to be 
positive definite, as is required by definition. The off-diagonal elements are not proportional 
to the current because each of the field that appears in the quadratic current is evaluated on 
a different solution. 

The conserved current also serves as an important check of our numerics. 

5.2.3 Numerical solutions and dispersion relations 

We could identify numerically three distinct quasi-normal modes in the mixed system. 

Figure (^) presents typical plots of the diagonal components of the Green's functions at 
vanishing k. Note that at zero k the equations of motion are decoupled and tx and et are 
identically zero. As a result the off-diagonal terms in the Green's function matrix is also zero. 

The first peak we found corresponds to the divergence in lm(Gxx) as u; — > 0. At fc = 
the equations of motion decouple and the divergence in lm{Gxx) corresponds precisely to the 
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Figure 7: The real and imaginary parts of the diagonal components of the Green's functions at 
d/C = 1.074, fc = 0. The dotted black curve corresponds to 1/w and the dotted red horizontal curve 
corresponds to 7r/2. (All plots are obtained at rn = 1-) 



dissipationless quasi-normal mode as computed in the previous section in the hydrodynamic 
limit. As already mentioned in section (^) however, this mode remains there with unaltered 
dispersion even in the finite d and q limit, where mixing becomes important. 

The second quasi-normal mode we found correspond to the second peak of lm{Gxx) in 
(0). Typical plots of it's dispersion are given in figure (P). The real part of the quasi-normal 
frequency begins roughly quadratic in k, but then asymptotically approaches lor = k. The 
imaginary part ajj decreases with increasing q, but then stabilises as k continues to increase, 
and show no sign of instabilities as far as we have checked. 

The first peak in Im(Gee) in fact corresponds to a pole at a; = — i(a -|- -Dfc^), for some 
constant a and D which depends on d and the embedding x- The peak gets shifted from uj = 
because it happens that ti approaches zero also at u = 0. This is the third quasi-normal 
mode we followed. A typical plot of the dispersion of this pole is given in figure (^). 

Beyond these three peaks visible in the plots of the spectral functions, they look com- 
pletely smooth and we have not attempted to identify more quasi-normal modes. 

Typical plots of the off-diagonal Green's functions at finite k are presented in figure ([lO|). 
It is not surprising that the imaginary part of Gxo goes negative since space {V) and time (T) 
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Figure 9: Typical plot of the dependence of uji on momentum q of the quasi-normal mode corre- 
sponding to the first peak in Gee- 



reversal symmetries are explicitly broken by the background electric field. This is consistent 
with the observation in [14] that positivity of ImGxe is guaranteed only when V and T are 
individually preserved. The diagonal components of the Green's functions at finite k look 
very similar to the k = case, except that the peaks are shifted, and therefore are not shown 
here. It is however important to point out that the imaginary parts of the diagonal elements 
remain positive definite at all frequencies. From the bulk perspective this is guaranteed by the 
conserved current, as already explained in ( 5.2.2| ). This is also required from the definition of 
the Green's function Go^o^(ti;) for any operator Oa in the dual field theory. The imaginary 
part of the Green's function is given schematically by 



IuiGoaOa (^) 
where Ea = Ea - /iiVa, E^fj = S 



'El 

a,f3 



<a\OA\P>\^ie 



±e-^='')5(u; + H„^) 



H^, Z is the partition function, and the sum is over all 
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Figure 10: The real and imaginary parts of tiic G^e at djC — 1.074, k = 0.7 and th = 1- 



eigenstates of the Hamiltonian. The + sign is taken for fermions and — for bosons. It is clear 
that every term in the sum is positive. 

We note that in the high frequency hmit lv ^ 1, the Green's function approaches that of 
the conformal result i.e. the vanishing baryon density d and temperature limit. The numerical 
result matches precisely the analytic result obtained in equation 

6. Fermionic fluctuations and spectral functions 

One of the most interesting aspects of a 1+1 dimensional system is the possibility of having 
spin-charge separation. Roughly speaking, this corresponds to charged fundamental fermions 
behaving as bound states of two degrees of freedom, one of which carries only spin and is 
called the spinon and the other which carries only charge and is called a holon. We use this 
as a motivation for studying the spectral function of world volume fermions. We however find 
no evidence of spin charge separation in these spectral functions. 

The analysis follows [51]. Firstly, in this setup, since we are considering a single D3 
probe, the fermions do not carry any U(l) charge. Hence the Dirac equation for fermions 
simply reads [52,53], 

(G + Ff%Z),V' = 0, (6.1) 

where 

D, = d, + ^LofTab . (6.2) 

Here Greek letters will indicate curved space indices while Latin indices will represent tangent 
space indices. 

For simplicity we will only consider the X = case leaving the massive case for future 
work. Following [51], the fermion spectral function exhibits two independent components 
given by Im^+ and Im^_ where ^± satisfy the following equations. 



= ±2DC± + i{A ±C) + i{A ± C)f± , 



(6.3) 




Figure 11: (a) A typical plot of the fermion spectral functions. The right branch in red is ^-|_ while 
the left branch in blue is (b) a vs djT for x = 0, r^f = 1. 



where 

„ . Frt^ ^ . 1 9rr , ^ ■ I 9rr , j-, drQxx /c a\ 

D = %—===, C = i. fc, A = i,- — - + Frt , , (6.4) 



^/ gxx\gtt\ V 9xx V \9tt\ '^gxx\fg rr\9tt I 

with 

Frt = ^=1= . (6.5) 

Here we have assumed a Fourier decomposition of the fermion fields with w,k denoting the 
energy and spatial momentum respectively. Infalling boundary conditions imposes 

e± = i , (6.6) 

at the horizon. It is easy to numerically solve for the spectral functions. A typical plot is 
shown in figure |Tl| (a). We should point out that there is no particle- hole asymmetry due 
to the absence of U(l) charge and also no sign of a fermi surface or a special momentum 
above which the quasiparticles disappear. The dispersion relation for large k is uj = zizk with 
decreasing widths. If we do a fit for the height of the peaks as a function of k of the form 
lmG{Q, k) = Ak", we find a trend as in figure |ll|(b) suggesting a significant deviation from 
the Landau Fermi liquid picture where a = —1 [51]. Curiously for large d/T, a seems to 
asymptote to a^, ~ 0.45. For low d/T, a is around 0.68. The exponent a/2 is the "anomalous 
dimension" of the fermion operators [15]. From this we can extract the parameter K for the 
Luttinger model which is given by 

a = ^iK + ^-2). (6.7) 

Using this and assuming a repulsive interaction in the effective Luttinger model (in effect 
choosing the root K < 1), we find K lies between 0.22 — 0.28. It is very interesting to 
compare the value for a for a massless fermion in the BTZ background, which can be extracted 
analytically. Considering the case that Tl = Tn and making use of the analytic Green's 
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function of a fermion in this background [54], and setting the fermion mass m = 0, we 
found that in the large k hmit the height of a QNM peak scales as k^R~^R = /c^/^, where 
= 3/4, hfi = 1/4 as defined in [54]. This gives a = 1/2 and K evaluates to K = 
2-V3^ 0.27, which is very close to what we found above numerically in the large d limit. 
Curiously, the value K = 0.28 is used to describe a wide range of experimental results in 
carbon nanotubes [55]\ 



7. Conductivities 



The conductivity can be extracted from the Green's function. At zero momentum, the con- 
ductivity is related to the current-current correlation function simply by 



a{w) 



w 



(7.1) 



whose temperature dependence at constant baryon density d can be readily obtained. Note 
that here we denote the dimensionful frequency by w, which is defined as u; = rH<jJ^^- The 
results are shown in figure (|l2|)^^. 




T 
C 



Figure 12: (a) and (b) show plots of the real and imaginary parts respectively of the conduc- 
tivities as a function of temperature at q = and lo = w/ru — 0.01/(7rT'), where d/C = 
1.04 (blue curve), 1.10 (black curve), and 1.20 (red curve). 

From figure (|l^), it is clear that the real part of the conductivity, becomes independent 
of T in the high temperature limit. Re {a) behaves like 



Re{a{w/d, W/C « l))|small T ~ Ke{a{w/d, w/C < l))|large T ~ T, 



(7.2) 



This can be compared with a normal Fermi liquid, whose real dc conductivity falls off as 1/T^ 
in the low temperature limit. It is also interesting to compare this with a simple version of 



^^The distinction between them is only important in this section where we would like to vary the temperature, 
instead of setting rn ~ 

^■^Note that we have not included the logorithmic dependence log rn in the imaginary part of the conductivity. 
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a spin half Luttinger liquid (see for example [15]), whose low frequency and low temperature 
conductivity can be parametrised in terms of the interaction strength, 

Re(adc)~T3-P''^p, (7.3) 

where 1/p is the band filling, and Kp characterises the interaction strength of the charge 
density waves [15]. 

Using the result for the Green's function in section (| 
frequencies, 

Re((T(7i;)) ~ tD^^, lm{a{w)) ^ . (7.4) 

Note that the \/w behaviour of lm.a{w) is deduced from the fact that Yle[Gxx) approaches a 
constant as u; ^ 0. This signals that the real part of the conductivity should contain a delta 
function 5{w) [14] on top of the l/nP' divergence, although the delta function would not show 
up in the numerics or a perturbative expansion in w, since w is never taken to be strictly 
zero. We again compare with the Luttinger liquid behaviour given by 

Re(cj(u;)) ~ w^^^"-^. (7.5) 

It is tempting to fit our results with those of the simple Luttinger liquid discussed above. 
At larger values of the temperature, where the conductivity approaches a constant, our results 
for the real part of the conductivity can be roughly fitted with the Luttinger liquid, with the 
parameter p^Kp ~ 3. Here Kp,p are treated as some effective parameters, a combined effect 
of all the degrees of freedom there in our theory. 

Typical plots of the conductivity vs frequency at constant temperature is given in figure 
(|T3|). Note that the real and imaginary parts are diverging at different powers of a; as u; — > 0. 



J5.2.1), we obtain also, at small 



(a) (b) 

Re[cr] Im[a-] 




Figure 13: (a) and (b) show typical plots of the real and imaginary parts respectively of the conduc- 
tivities as a function of the frequency aX q — ru — ^■ 
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8. Discussion 



We studied 1+1 dimensional holographic quantum liquids using a D3-D3 defect CFT. We 
set up the thermodynamics and explained how the boundary terms differ when compared to 
higher dimensions. Unlike higher dimensions, it is necessary to introduce a gauge potential 
counterterm to remove the log divergence arising from it. We explained how to compute 
both bosonic and fermionic spectral functions in this setup. Unlike higher dimensions we 
found a dissipationless mode and a purely dissipative mode. The dissipationless mode is a 
pole in the density fluctuations, and can thus be interpreted as a sound mode. This is what 
is generally found in a Luttinger model after diagonalisation by bosonisation, from which one 
obtains the charge density wave as an eigen-fluctuation of the Hamiltonian, which disperses 
at a particular velocity v depending on the coupling of the mother fermions [15,21]. In our 
case, the velocity is simply v = 1, taking the conformal value^^. 

To the best of our knowledge, the latter purely dissipative mode is not present in any 
Luttinger model. It could be that the appearance of this mode is a large N artifact. This 
issue clearly merits further study. In the extended Hubbard model in 3+1 dimensions [56], a 
purely dissipative spin mode and a non-dissipative charge mode was found. We are not aware 
of similar studies for the extended Hubbard models in 1+1 D. In our case, when d = 0, the 
dissipative mode disappears and it is likely that a Luttinger description holds and one may 
be able to extract the parameters v and K used in this description. We find that v = 1 in 
all cases. For the fermions, assuming an effective Luttinger description, it was found that 
the parameter K > 0.22. In fact for a wide range of d/T we found K ~ 0.28. This value is 
reported in the literature [55] as describing a wide range of phenomena in carbon nanotubes. 
A simple estimate of K from the Luttinger model, making explicit use of the graphite band 
structure gives K ~ 0.2 [57], which demonstrates that K is a material dependent quantity. In 
fact more recently studies of Gallium Arsenide quantum wires gives a value of K ~ 0.66 [58]. 
These results, combining with the plot of a vs d/T in figure (11) would suggest that there 
is an upper bound on a, and thus a lower bound^"^ on K > 0.22. Clearly this deserves a 
closer look. It should be possible to determine the effective K and v for the bosonic sectors 
in the probe brane model but for this one needs to work out the fermionisation map between 
the Luttinger fermions and the AdS bosons. Moving the branes away from the origin of 
the moduli space [25] may also lead to different physics. A model that will be interesting 
to study is the D3-D7 setup making a 1+1 dCFT. This model will have chiral fermions. It 
would be interesting to compare with what is known about chiral Luttinger models [64]. A 
first step in this direction has been taken in [65]. Moreover, it will be interesting to answer: 
What is the holographic analogue of spin charge separation? In the spin-half Luttinger liquid, 
the charge and spin degrees of freedom are identified with n{k)^ ± n{k)^ respectively, where 

^•^This is probably a result of working in the probe limit, where we are breaking conformal invariance 
"weakly" . Deviation from the conformal result is probable upon taking into account back reaction. 

^■'One of the achievements of AdS/CFT is to come up with the famous viscosity bound 77/s > l/47r [59-61]. 
Although the precise number is doubtful [62,63], it is thought that a bound exists. In a similar spirit, one can 
look for other physical quantities where AdS/CFT methods would lead to a bound. 
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n{k)i^l correspond to the Fourier components of the density operator of the spin up and down 
fermions respectively [15]. We have studied in this paper the Green's function of the charge 
density of the world- volume U{1), which we can naturally identify as the Green's function of 
the charge density wave, whose dual operators involve the bilinears of fermions belonging to 
the hypermultiplets [25]. We will treat the R-symmetries indices of the fermions in the dual 
dCFT as spin labels. One could then look for suitable fluctuations on the probe brane which 
would allow us to compute the Green's function of the spin-density wave. A candidate of 
such a fluctuation is that of A,^-^ and one could generally expect that the quasi-normal modes 
would disperse differently from the charge density degree of freedom [66], leading to a natural 
realisation of spin-charge separation. So far we have identified the lowest mode of A^^-^ which 
is also a purely dissipative mode. This is already suggestive of a spin-charge separation, with 
the spin density wave dispersing at v = and the charge density wave at v = 1 [66]. 

One feature that exists in Luttinger liquids is the Schottky anomaly in the specific heat 
[67]. This happens when the number of states available in the microcanonical ensemble is 
finite and results in a bump in the specific heat. In Luttinger liquids, this is expected to 
happen at low temperatures. Although no such bump happens for the x = case, we have 
not been able to convince ourselves about the accuracy in our numerics to rule this out at 
X 7^ 0. Moreover, as shown in appendix A, for small Xi d one can have an analytic handle on 
the specific heat. We find there that for non-zero Xi the difference in the specific heat from 
the zero x result is proportional to cP and is positive. Since for large temperatures, both cases 
will lead to ~ irT, it seems that one prerequisite to have a bump is that the change in the 
specific heat from the x = case is positive! So it could still be that a better handle on the 
numerics will lead to something interesting like the Schottky type specific heat. 

g^cv We would like to conclude with the following 

observation. Since the low temperature specific 
heat behaves like Cy oc and the high temper- 
ature specific heat behaves like oc T, it is clear 
that in the intermediate regime, dxCv should have 
an interesting feature. As shown in figure (|l4|), we 
indeed find that for x = Oi dxc^ exhibits a max- 
imum for T = \f2dl'K. As we turn on x ex- 
pect the position of this maximum to get shifted to 
Figure 14: drCv vs T jd for x = 0. The the left (since generically the combination \fSF—(? 
maximum occurs at T = \f2dl tt. appears in place of d) and eventually vanish near 

d ^ c. The fact that the low temperature specific 
heat behaves like and the high temperature one behaves like with h < £ appears quite 
generic in these models. Hence we expect this maximum in dxCy to persist in all such defect 
CFTs. 

After this work was completed [68] appeared which deals with related issues. 
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A. Thermodynamics at small d and x 

We have discussed the thermodynamics of the system at both vanishing and finite x- We find 
that at finite x numerical solutions of the embedding are used in the computation. We can in 
fact extract further analytic result in the limit of small d and x, and expand the bulk action 
up to quadratic order in d and x- 
Doing that, we have 

Er = -, Air) = dlog—, (A.l) 
r th 

which means /i = /I + d\og rn = 0. 
The solution for x is given by 

, , C^[r(l/4)^Fi(l/4,l/4,l/2,n^) -2t/^r(3/4)^Fi(3/4,3/4,3/2,^^)] 

4rj,r(l/4)r(3/4) ' 

where we have fixed the boundary condition at the horizon such that it is regular, and that 
at the AdS boundary the coefficient of the leading log term is simply c = C/rn- 

We can substitute these expressions into the action and the counter terms. This gives 

— = rtri logr/7 , 

iVo 2 4 2 ^ 

S , (d^-c^), , {d^-(?), ^ 



which agrees with the small d expansion of the analytic result ( 3.15 , 3.14 ) at zero c. 

The expression (d^ — (?) is strictly positive in the stable region and negative in the 
unstable region. 



B. The mixed modes Lagrangian 

We would like to specify in the table below the leading near horizon and boundary behaviour 
of the various functions we defined in equation (|5.1| ). 
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u — > 00 
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dxo^/l-Xo 
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cd log u 
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Si 
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U 
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S3 
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(B.l) 
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C. More on stability 



Singularities of the spectral functions generally occur in the complex w plane. When a; 
is complex, 0_fc and similarly a^-k cease to be the complex conjugates of Qk and a^^k 
respectively. We can, however, still obtain some analytic results on the stability of the system 
(i.e. determine the sign of uji at which the singularity is located). To that end, let us take 
equation ( |5.4| )multiplied by at^^i and subtract it by ax^k multiplied by the complex conjugate 
of ( [5. 3D and finally adding to them Ok times ( ^.12 ). The resultant equation is given by (we 
will omit k in the expression) 

2atdu{S2a[) - laMSsa'x) + 2ea„(r4e') + (T^ - 2(ri + k^T2 + u;^n))\Q\^ 
—2k'^ Si\at\'^ + 2u)'^ Siluxl"^ + 2kdtaxSi{u! — uj) 
+Dik{ojQax - katQ) - atdu{D^Q + D2Q') 

-Dik{kaS + + ©(^30! - 5(^2al)) = 0, (C.l) 

where the primes also denotes partial derivatives w.r.t. u. Taking the imaginary part and 
integrating by parts, we are left with 

\m[d{2atS2a[ - 2axS^a'x + 2er4e' + atD^ + ai^2e' + BZ^saO] 

= -a;/(4L^/j,(r3|e|2 + Si\ax\'^) + Dik{axQ + a^O) + 2kSi{atax + a^a*)), (C.2) 



where uj = ujr — iu>j. The cross terms on the r.h.s. of (|C.2| )can be replaced by a single term 
using (|5.3| ). After some simplification we are finally left with 

du{Im(2atS2a[ - 2a^53a' + 2eT4e') + ^Re(2wa^53a' )} 

\^\ 

= -u;i{Au;rTs\Q\^ - 253|<P^). (C.3) 

1^1 

We have dropped cross terms inside the partial derivatives since we will eventually inte- 
grate over u from the horizon to the boundary, and the D terms vanish at both limits. The 



expression can be further simplified using ( 5^ , 0| ), and explicitly integrating over u from the 
horizon to the boundary, we have 



I a; - A; J 

00 



= uJiujR V{u), (C.4) 

where 

V{u) = A{-Ts)\ef + 2sJ-^. (C.5) 

\^\ 

We have also made used of the explicit form of the solutions of Ex and G at the horizon. 
From the explicit form of T3 which is negative definite and which is positive definite, we 
can conclude that V{u) is positive definite. 
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Since the quasi-normal modes appear at the same time in Q and E^, the first term 
evaluated at the boundary is zero at a quasi-frequency. The l.h.s is evaluated explicitly at 
the horizon, and is positive definite. Therefore it implies that loj has to be positive, at the 
level of our approximations, where these perturbations are small. The analysis here does not 
guarantee stability, since the same analysis would have been valid in the case of a D7 probe, 
where we are well aware of the presence of instability from [28,36,37]. It is important to 
note that if the quasi- normal mode occurs at ujr = 0, the argument breaks down and our 
equations of motion could produce a quasi-frequency where ojj takes a negative value. In 
fact we found precisely such a mode which gives a negative ujj where the heat capacity turns 
negative, indicating instability. 

D. The explicit form of the Green's functions from the probe D3-brane 

We would like to work out explicitly the form of the Green's functions, including the counter 
terms for the probe D3-branes along the lines outlined in section ( |5.2.1| ). 

Consider the probe D3 in the AdS black hole background. Ignoring for the moment the 
correct normalization, let us write the solutions in the boundary limit as 

= ei log[n] + eo + O(^), = ^ + O(^), 
i. = to + 0(^),t. = t.^ + ^ 

(D.l) 

and recall that 

= Eoct + Tott, Ex = EqCx + Totx 

For notational simplicity, we denote {Q-k, Ex^k} — {^k: Exk}- This relation however, holds 
literally only when to is real. We have 

which explicitly evaluates to 

Mee = |e/plog[noo] +eoe/ - |eip. 
Met = - (eiiU - ti) - eiU log[noo] - {Uei log[uoo] + hei) + toe/) , 

Mte = Wt, 

M-tt = -\ti log[noo] + tip + logfuoo] + htu (D-3) 

so that the Green's function matrix again appears to be divergent. Here we have absorbed a 
factor of — k'^) into ex and tx momentarily to avoid clumsy notation. 
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At the same time, counter terms have to be taken into account. Let us therefore evaluate 
the counter terms in (|3^,3.5), where the quadratic terms are given by 



C = / (fx{ul^QQ{l 



1 



(D.4) 



Using the equations ( p^j5.8|) , and the asymptotic behaviour of the functions Si^2 ~ ±n/2 
and that the Di terms are subleading at the boundary, we have, 

ujEx(u,k) Q —kEx(u,k) n 
ax{u ^ oo,k) = — -^ + ax, at[u ^ oo, k) = —-^^ j-^ \-a^, (D-5) 



2 



for some constants a?.. 



Combining the above equations we have 

kat{u — > oo. A;) + u!ax{u — > oo. A;) = Ex{u — > oo, k) + uja^ + ka^. (D.6) 



This imphes that 



a;a° + ka"! = 0. 







(D.7) 



As a result, at the boundary, the counter terms A^A^ evaluates to ExEx/{oj'^ — k"^)- 
The contribution of the counter terms to the Green's functions are thus given by 



|e/p(log[noo] + logr/^) - {e^ei + eoe/) + |eip. 



1 



Get = ^ ((i/ei log[uoo] +tiei) - t^ei - t/ei - ei{ti - ti) - eito + eit; log[uoo]) , 

Cte = Cet, 

Clt = \tl\og[Uoo\+tl\'^ - {\ti\'^{\og[uoo\+\ogrH)+titi+titi). 



(D. 



Note that the logrn contribution is explicitly real, which we will drop from now (.i.e. setting 

TH = 1). 

Combining with (p.2|), and restoring factors of — k'^ we have 



Gee Get 
Gte Gft 



( 26; eg 



Putting back the correct normalisations as defined in ( 5.17 , 5.18 ), we have finally 



Gee Get 
Gte Gtt 



eo 



to 



ei 



to 



ti{uj-^-k-^) ^ ei 



^ti 



as advertised in ( 5.22| ). 
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